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Abstract 

Using an eikonal analysis, we simultaneously fit a QCD-inspired parameterization of all accelerator 
data on forward proton-proton and antiproton-proton scattering amplitudes, together with cosmic ray 
data (using Glauber theory), to predict proton-air and proton-proton cross sections at energies near y/s w 
30 TeV. The p-air cosmic ray measurements greatly reduce the errors in the high energy proton-proton 
and proton-air cross section predictions — in turn, greatly reducing the errors in the fit parameters. From 
this analysis, we can then compute the survival probability of rapidity gaps in high energy pp and pp 
collisions, with high accuracy in a quasi model-free environment. Using an additive quark model and 
vector meson dominance, we note that that the survival probabilities are identical, at the same energy, 
for 7p and 77 collisions, as well as for nucleon-nucleon collisions. Significantly, our analysis finds large 
values for gap survival probabilities, ~ 30% at i/s — 200 GeV, ~ 21% at \/s = 1.8 TeV and m 13% at 
Js = 14 TeV. 
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Rapidity gaps are an important tool in new-signature physics for ultra-high energy pp collisions. In this 
note, we will use an eikonal model to make a reliable calculation of the survival probability of rapidity gaps 
in nucleon-nucleon collisions. 

In an eikonal model^], we define our (complex) eikonal x(b, s) so that a(b,s), the (complex) scattering 
amplitude in impact parameter space 6, is given by 

a ( bj s ) = i (i _ e *x(M)) = * fi _ e -x,(M)-Mx R (M)) (1) 
Using the optical theorem, the total cross section atot(s) is given by 

0tot(s) = *j [l ~ e-^( fe ' s ) cos(x a (&,«))] d 2 b, (2) 
the elastic scattering cross section <7 e i(s) is given by 

TclasticOO = / |l-e-Xl< 6 >">-*XR< 6 '">| 2 d 2 b (3) 

and the inelastic cross section, <Ti ne iastic(s), is given by 

^inelastic (s) = - . elastic ( S ) = / [l - e^] d%. (4) 

The ratio of the real to the imaginary part of the forward nuclear scattering amplitude, p, is given by 

Re (i( f 1 - e-XrM+ix^,*)) 

P(s) = > Jr (5) 

Im - e-Xj^.^+ixR^s)) jpU 

and the nuclear slope parameter B is given by 

B = jb 2 a(b,s)d 2 b 
2 J a{b,s)d 2 b' 

It is readily shown, from unitarity and eq. (||), that the differential probability in impact parameter space 
b, for not having an inelastic interaction, is given by 

p . , .. = p-2x r (M) (7) 
A no inelastic o . yij 

Our even QCD-Inspired eikonal Xevcn is given by the sum of three contributions, gluon-gluon, quark-gluon 
and quark-quark, which are individually factorizable into a product of a cross section cr(s) times an impact 
parameter space distribution function W(b;p), i.e.,: 

X cvcn (s,b) = Xgg(s,b) + Xqg(s,b) + Xqq {s,b) 

= i [<Tgg(s)W(b ; ^ gg ) + a qg (s)W(b ; v / ^W%g) + o- qq (s)W(b ; p qq )] , (8) 

where the impact parameter space distribution function is the convolution of a pair of dipole form factors: 

W{b-^) = -£ r {»bfK^b). (9) 
967T 

It is normalized so that J W(b ; p)d 2 b = 1. Hence, the cr's in eq. (||) have the dimensions of a cross section. 
The factor i is inserted in eq. (||) since the high energy eikonal is largely imaginary (the p value for nucleon- 
nucleon scattering is rather small). The total even contribution is not yet analytic. For large s, the even 
amplitude in eq. (||) is made analytic by the substitution (see the table on p. 580 of reference ||, along with 
reference Q) s — > se~ l7T / 2 . The quark contribution x qq (s,6) accounts for the constant cross section and a 
Regge descending component (oc 1/s/s), whereas the mixed quark-gluon term Xqg( s i b) simulates diffraction 
(oc logs). The gluon-gluon term x gg (s,6), which eventually rises as a power law s e , accounts for the rising 
cross section and dominates at the highest energies. In eq. (g), the inverse sizes (in impact parameter space) 
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/x qq and /i gg are to be fit by experiment, whereas the quark-gluon inverse size is taken as ^//i qq /i gg . For more 
detail, see ref. (2j. 

The high energy analytic odd amplitude (for its structure in s, see eq. (5.5b) of reference M, with a = 0.5) 
that fits the data is given by 

X° dd (M = -<?odd W(b lf i odd ), (10) 

with cr dd oc l/v^j an d with 



W(b, Modd) = ^(Moddfo) 3 ^ 3 (/iodd6), (11) 

normalized so that J W(b ; fJ- dd)d 2 b — 1. Hence, the tr dd in cq. ( |Io| ) has the dimensions of a cross section. 
Finally 

Xpp — Xcvcn i Xodd- (12) 

We have constructed a QCD-inspired parameterization of the forward proton-proton and proton-antiproton 
scattering amplitudes pi which is analytic, unitary, satisfies crossing symmetry and, using a x 2 procedure, 
fits all accelerator data|6| of ertot, nuclear slope parameter B and p, the ratio of the real-to- imaginary part 
of the forward scattering amplitude for both pp and pp collisions; see Fig.|l| and Fig.^j. In addition, the high 
energy cosmic ray cross sections of Fly's Eye [Q and AGASA Q experiments are also simultaneously usedj|]. 
Because our parameterization is both unitary and analytic, its high energy predictions are effectively model- 
independent, if you require that the proton is asymptotically a black disk. Using vector meson dominance 
and the additive quark model, we find further support for our QCD fit — it accommodates a wealth of data 
on photon-proton and photon-photon interactions without the introduction of new parameters^. A major 
difference between the present result, in which we simultaneously fit the cosmic ray and accelerator data, 
and our earlier result in which only accelerator data are used, is a significant reduction (about a factor 
of 2.5) in the « 1.5 % error of our prediction for a pp at = 30 TeV, which results from major reductions 
in our fit parameters. 

The plot of (Tpp vs. y/s, including the cosmic ray data that have been converted from cr^air 1° a vv^ ^ s 
shown in Fig. [l]. Clearly, we have an excellent fit, with good agreement between AGASA and Fly's Eye. The 
overall agreement between the accelerator and the cosmic ray pp cross sections with our QCD-inspired fit, 
as shown in Fig.|l], is striking. 

As an example of a large rapidity gap process, we consider the production cross section for Higgs-boson 
production through W fusion. The inclusive differential cross section in impact parameter space b is given by 
Js= = ctww^h W(b ; (iqq), where we have assumed that W(b ; /i qq ) (the differential impact parameter space 
quark distribution in the proton) is the same as that of the W bosons. 

The cross section for producing the Higgs boson and having a large rapidity gap (no secondary particles) 
is given by 

= ^h^iM^'^" = ^ww^h (13) 
d 2 b db 2 



In eq. (13) we have used eq. (Q) to get the exponential suppression factor, and have used the normalized 

impact parameter space distribution W(b;fi qq ) = ^(/i qq 6) 3 K 3 ()j, qq b), with fi qq = 0.901±0.005 GeV, whose 
error comes from the fitting routine)^]. 

We now generalize and define < IS*! 2 >, the survival probability of any large rapidity gap, as 

< \S\ 2 >= [ W(b ; n qq )e~ 2 ^ s ' b U 2 b. (14) 



We note that the energy dependence of the survival probability < IS*! 2 > is through the energy dependence 
of Xn the imaginary part of the eikonal. 

For illustration, we show in Fig. |3|a plot of Im%p P and the exponential damping factor of eq. (O), as a 
function of the impact parameter 6, at -Js — 1.8 TeV. The results of numerical integration of eq. (|14[) for the 
survival probability < S 2 > at various c.m.s. energies are summarized in Table |l|. As emphasized earlier, the 
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C.M.S. Energy (GeV) 


Survival Probability (pp), in % 


Survival Probability {pp), in % 


63 


37.0 ±0.9 


37.5 ±0.9 


546 


26.7±0.5 


26.8 ±0.5 


630 


26.0 ±0.5 


26.0 ±0.5 


1800 


20.8 ±0.3 


20.8 ±0.3 


14000 


12.6 ±0.06 


12.6 ±0.06 


40000 


9.7 ±0.07 


9.7 ±0.07 



Table 1: The survival probability, < |>S| 2 >, in %, for pp and pp collisions, as a function of c.m.s. energy. 



errors in < l^l 2 > are quite small, due to the fitting parameter errors being very small when we determine 
the eikonal using both accelerator and cosmic ray data. 

Further, we find for the quark component that the mean squared radius of the quarks in the nucleons, 
< i? 2 n >, is given by < R^ n >= j b 2 W{b; ^ qq ) d 2 b = 16/^ qq = 19.70 GeV -2 . Thus, 6 rms , the r.m.s. impact 
parameter radius is given by & rms = 4//i qq = 4.44 GeV -1 . Inspection of Fig. || (1.8 TeV) at 6 rms shows a 
sizeable probability for no interaction (e~ 2x i) at that typical impact parameter value. 

In ref. 0, using the additive quark model, we show that the eikonal x yp f° r IP reactions is found by 



substituting a 
substituting a 

x 7P M) = * 

and 



§CT, fl 
§CT, (1 



2 H into x eve "( s i b), given by eq. 
|/x into x lp (s,b). Thus, 



|). In turn, x~ n f° r 77 reactions is found by 



2 /3 2 /3 2 /3 

-a ge (s)W{b; y ^Mgg) + ^qs( s ) W ( b ; d -^MqqMgg + ^o- qq {s)W (b 5 V 



4 3 4 3 4 3 

~<T gg (s)W(b; -/x gg ) + gCr qg (s)PF(6; -VMqqMgg + gCr qq (s)PF(6; -/x qq ) 



X 77 («,6) = 

It can be shown jl[], that by an appropriate change of variables, that 



<|^| 2 >= j W{b-^-^e-^( P ^dH 
where x 7P (s, b) is given by eq. (|l5|), and 



IF(&;^ qq ) e —r^^ 2 6, 



< |S 77 | 2 >= 

where x 77 ( s i 6) is given by eq. (|16[), are 6oi/i equal to 

< |S m f >= J PF(6 ;Alqq )e" 2x r"(^ b )d 2 6 

where x cven ( s i^) is given by eq. (||). Thus, we have the interesting result that 

< \s™\ 2 >=< is" 77 ! 2 >=< is 11111 ! 2 > . 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



Neglecting the small differences at low energy between pp and pp collisions, we see from eq. ( p0| ) that 
< \S\ 2 >, the survival probability for either nucleon-nucleon, jp and 77 collisions, is reaction-independent, 
depending only on </s, the cms energy of the collision. We show in ref. |Q that this result is true for any 
factorization scheme where the eikonal factorizes into sums of (7i(s) x Wi(b; /i), with the scaling feature that 
the product cr^/x 2 is reaction-independent — not only for the additive quark model that we have employed 
here. The energy dependence of the large rapidity gap survival probability < l^ 2 > calculated from eq. ( |l9| ) 
is given in Fig. 0. 
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This somewhat surprising result can be more readily understood once one realizes that the survival 
probability is a function of the product a/j, 2 and the dimensionless variable x — /i qq 6. This is most simply 
seen in a toy model in which the even eikonal is given by 

even, ,\ ttz/u \ 2 (Vqq b ) 3 K 3(Vqqb) 2 X 3 K 3 (x) 

V° von (s, 6)toy = OqqW{b- Mqq ) = a qqt i 2 qq x 96 ^ = <W qq x Q6lT . ( 21 ) 

and hence, 

< l-SHtoy >=7^J e~ 2a ^<^ K3{x)/9 ^d 2 x. (22) 

In our scheme OqqP 2 ^ is the same for nucleon-nucleon, jp and 77 interactions, since cross sections are 
nd I for 7p and 77, respectively, whereas fj, 2 ^ is divided by ^/| and |. Although x is 

different for the 3 processes — being fi q b for nucleon-nucleon, \J\p-qqb f° r IP an d f/^qq^ f° r 77 — it only plays 
the role of an integration variable and therefore the dependence on the subprocess label disappears. We see 
from eq. (|2^) that S=S(CT qq /i? q ), which is process-independent. This argument is easily generalized to the 
full eikonals of eq. (||), eq. ( ]l5[ ) and eq. (pf). 

The physics is now clear. One could be tempted to conclude that the survival probability is larger for 
7P than for pp interactions because there are only 2 quarks in the photon and 3 in the proton to produce 
additional inelastic collisions filling the gap. And that is true and it is reflected by the factor 2/3 change 
in the cross section. This is not the whole story, however. In the eikonal model, the inverse transverse 
size of the 2-quark system (photon) is larger than that of the proton (the |/i qq factor) and the two effects 
compensate. Therefore, we have the same < IS*! 2 > for all processes. 

We have been able to calculate the survival probability < \S\ 2 > to a high degree of accuracy by using an 
eikonal that has been found by fitting both accelerator and cosmic ray data over a very large energy range 
which includes the LHC energy. Our numerical results are considerably larger than other calculations ||lo|, O , 
|l2| . In the case of ref . Jl(| and ref. [E2I it is probably due to their using a Gaussian probability distribution 

in impact parameter space, whereas our distribution, W(b,p qq ) — g^(n qq b) 3 K 3 (fj, qq b) which is the Fourier 
transform of the square of a dipole distribution, has a long exponential tail e~' tb , significantly increasing 
the probability of survival. In the case of ref. the authors determine the parameters for their minijet 
model using only the Tevatron results. Our large values are more in line with the earlier predictions of 
Gotsman et al^j for what they called Regge and Lipatovl and Lipatov2 models, although with somewhat 
different energy dependences than that shown in Table 1. The color evaporation model of Eboli et aZ[[l4} gives 
somewhat larger values than our calculation, but again with a different energy dependence. Most recently, 
Khoze et al^jfy, using a two-channel eikonal, have calculated the survival probabilities for rapidity gaps in 
single, central and double diffractive processes at several energies, as a function of b, the slope of the Pomeron- 
proton vertex. For double diffraction, they have a large range of possible parameters. Choosing 2b = 5.5 
GeV -2 (corresponding to the slope of the electromagnetic proton form factor), they obtain < \S\ 2 >— 0.26, 
0.21 and 0.15 at y/s = 0.54, 1.8 and 14 TeV, respectively. These survival probabilities are in excellent 
agreement with our values given in Table 1. However, their calculations for other choices of 2b and for single 
and central diffractive processes do not agree with ours, being model-dependent, with their results varying 
considerably with their choice of parameters and model. 

We see that there is a serious model dependence, both in the size of the survival probabilities and in 
their energy dependence. Further, until now, there has been no estimates for gap survival probabilities for 
7P and 77 reactions. Thus, we hope that our quasi model-independent fit to experimental data on pp and pp 
total cross sections, p values and nuclear slopes B, over a large energy range, y/s = 15 GeV to 30,000 GeV, 
provides a reliable quantitative estimate of the survival probability < \S\ 2 > as a function of energy, for 
both pp, pp 7 7p and 77 collisions. The fact that our estimates of large rapidity gap survival probabilities are 
independent of reaction, thus being equal for nucleon-nucleon, 7p and 77 processes — the equality surviving 
any particular factorization scheme — has many interesting experimental consequences. 
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Fi gur6 1: The fitted app and <Tp P , in mb vs. y/s, in GeV, for the QCD-inspired fit of total cross section, B and p for both pp 
and pp. The accelerator data (squares are pp and circles are pp ) and the cosmic ray points (diamond, Fly's Eye and triangles, 
AGASA) have been fitted simultaneously. The cosmic ray data that are shown have been converted from cri 1 ^^ to a vv . 




Figure 2: The fitted values for the nuclear slope parameters B pp and B pp , in (GeV/c) 2 vs. y/s, in GeV, for the QCD-inspired 
fit are shown in (a). In (b), the fitted values for pp p and p vp are shown. 
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Figure 3: A plot of the eikonal Im\ and the exponential damping factor e 2 * Im x f or pp collisions, at i/s = 
impact parameter b, in (GeV/c) . 
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Figure 4: The energy dependence of < |S| 2 >, the large rapidity gap survival probability vs. y/s, in GeV. 
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